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Figure 1. Iso-mass and iso-composition lines for neutralino in the M^-pt plane for tan [3 = 8. 
Dashed lines are lines of constant m x ; solid lines are lines of constant gaugino fractional weight P. 
The upper figure refers to the case jj, > 0, the lower to jj, < 0. 
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Figure 2. Neutralino relic abundance £! x /i 2 as a function of m x for the three representative 
neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed line), 0.9 (long-dashed line). 
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Figure 3. Same as in Fig. 2, except that now composition are: P = 0.01 (dotted line), 0.5 
(short-dashed line), 0.99 (long-dashed line). 
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Figure 4. Scatter plot for the neutralino relic abundance £i x h 2 as a function of m 
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Figure 5. tanh 2 (t / ta) and Ta are given as functions of m x for positive values of the parameter 
H, for the three representative neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed line), 
0.9 (long-dashed line). This figure refers to the Earth. 
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Figure 6. Same as in Fig. 5 for the case of the Sun. To make the comparison of the annihilation 
rates for the Earth and for the Sun easier, the plot shown here for the Sun refers to the effective 
annihilation rate, defined as Ta times the ratio (di/d'j) 2 = 1.8 x 10 -9 , where d\ is the radius of the 
Earth and c?2 is the Sun-Earth distance. 
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Figure 7. Fluxes of the up-going muons as functions of m x for x~X annihilation in the Earth, 
for the three representative neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed line), 
0.9 (long-dashed line). The threshold for the muon energy is = 2 GeV. 
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Figure 8. Same as in Fig. 7 except now the compositions are: P = 0.01 (dotted line), 0.5 
(short-dashed line), 0.99 (long-dashed line). 
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Figure 9. In the M^-pt plane are shown the regions explored by Kamiokande (full circles) and 
the domains explorable by an improvement factor of 10 (asterisks) and of 100 (dots) in sensitivity. 
Dashed lines denote iso-mass curves with m x = 30, 100, 300, 1000 GeV. Solid lines denote curves of 
iso-composition (here P = 0.01, 0.1, 0.5, 0.9, 0.99 are shown). 
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Figure 11. Event rates for detection of neutralino dark matter for a (natural composition) 
Germanium detector for the three representative neutralino compositions P = 0.1 (dotted line), 0.5 
(short-dashed line), 0.9 (long-dashed line). The rate refers to the events with an electron-equivalent 
energy in the range (2-4) KeV. 
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Figure 12. Same as in Fig. 11 except that now the three representative neutralino compositions 
are P = 0.01 (dotted line), 0.5 (short-dashed line), 0.99 (long-dashed line). 
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Figure 13. Explorable regions in the M^-pt plane for neutralino-Ge interaction. Asterisks 
denote regions where the integrated rate i?; nt is in the range 0.1i? e xp < -Rint < Rexp, Rexp being 
the present experimental bound. Dots denote regions where the integrated rate i?; nt is in the range 

0.01i?exp < -Rint < 0.1-Rexp, 
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ABSTRACT 



We review some properties of the neutralino as a candidate for dark matter in 
the Universe. After presentation of evaluations for the neutralino relic abundance, 
possibilities for its direct and indirect detections are discussed, with emphasis for 
measurements at neutrino telescopes. 



1 The cold side of the dark: neutralinos 



As is widely discussed in the literature, neutralino may be considered as a 
very natural candidate for Cold Dark Matter (CDM). This property rests on two 
assumptions, i.e. that: i) R-parity is conserved, ii) neutralino is the lightest su- 
persymmetric particle (LSP). Very convenient theoretical frameworks where dark 
matter neutralino phenomenology may be easily studied are provided by the Min- 
imal Supersymmetric Standard Model (MSSM) and by its implementation in a 
Supergravity theory (SUGRA) [1-2]. 

The neutralino (%) is defined as the lowest-mass linear combination of 
photino, zino and higgsinos 

X = a i7 + a 2Z + a 3 H® + a±Hl . (1) 
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Here 7 and Z are linear combinations of the U(l) and SU(2) neutral gauginos, B 
and W3, 

7 = cos 9wB + sin OwWs, 
Z = — sin 9wB + cos OwWs, 

and 0w is the Weinberg angle. 

The neutralino mass m x and the coefficients a,{ depend on the parameters: 

fj, (Higgs mixing parameter), Mi, M2 (masses of B and of W3, respectively) and 
tan/3 = v u /vd (v u and are the v.e.v.'s which give masses to up-type and down- 
type quarks). It is customary to employ the standard GUT relationship between 
Mi and M 2 : Mi = (5/3) tan 2 W M 2 ~ 0.5 M 2 . We use this assumption here. 

In the following for the parameters M 2 and fj, we will consider the ranges: 20 
GeV < M 2 < 6 TeV, 20 GeV < \fx\ < 3 TeV. tan/3 will be taken at the representative 
value tan/3 = 8. 

For the evaluation of the neutralino relic abundance and of the event rates for 
direct and indirect neutralino detections one also has to assign values to the masses 
of a large number of particles, namely to the Higgs bosons and to the Susy scalar 

partners of leptons and quarks: sleptons (I) and squarks (q). In the MSSM scheme 
these values are assigned arbitrarily: a standard procedure consists in assuming 
mass degeneracy both for sleptons and for squarks. Only for stop particles of 
different helicities non-degeneracy is sometimes introduced, since, under certain 
circumstances, stop-mixing may generate sizeable effects [3]. As for the neutral 
Higgs bosons we recall that in the MSSM there are three neutral Higgs particles: 
two CP-even bosons h and H (of masses nth, ran with ran > rrih) and a CP-odd 
one A (of mass tua)- Once a value for one of these masses (say, rrih) is assigned, 
the other two masses (m^, ran) are derived through mass relationships depending 
on radiative effects. 

Implementation of MSSM with supergravity sets a much more constrained 
phenomenological framework, since SUGRA establishes strict relations between all 
the masses in play and the few fundamental theoretical parameters: A and B (ap- 
pearing in the soft symmetry-breaking interaction terms), mo (common scalar mass 
at the GUT scale), m 1 / 2 (common gaugino mass at the GUT scale) and fj,. Further- 
more, other specific theoretical requirements (features of the symmetry breaking, 
condition that the neutralino be the LSP, ...) strongly restrict the whole parameter 
space. This has important consequences; for instance, it constrains the neutralino 
to compositions with dominance of the gaugino components. 

In the present note we wish to discuss the event rates for neutralino searches 
in a framework which is not too much constrained by theoretical requirements; we 
then simply adopt here a MSSM, with choices for the free parameters that are only 
restricted by experimental bounds. Our main concern is to discuss the minimal 
sensitivity required in experimental devices in order to undertake a significant in- 
vestigation of neutralino dark matter. For this reason we present here evaluations 
where the unknown masses are assigned the smallest values compatible with exper- 
imental lower bounds; this usually provides maximal values for the signals. To be 



(2) 



definite, in the following we will set the sfermion masses at the value mj = 1.2 ra x , 

when m x > 45 GeV, mj = 45 GeV otherwise. Only the mass of the stop quark is 

assigned a larger value of 3 TeV. Higgs mass rrth is set at the value of 50 GeV. The 
top mass has been fixed at m t = 150 GeV. 

We show in the M2 — p> plots of Fig.l the iso-mass lines and the iso- 
composition lines for the neutralino. Along an iso-composition line the compo- 
sition parameter P, defined as the gaugino fractional weight, i.e. P = a\ 2 + a2 2 , 
is kept fixed. Shown in Fig.l are the iso-composition lines referring to the values 
P = 0.01,0.1,0.5,0.9,0.99 (i.e., from a very pure higgsino composition to a very 
pure gaugino composition). 

In the following we will characterize a neutralino state by the values of m x 
and P; we will accordingly present our results for the relic abundance and the event 
rates in terms of these two parameters. 

2 How many neutralinos around us? 



For the computation of the direct and indirect event rates for neutralino 
one has to use a specific value for the neutralino density p x . Obviously, it would 
be inappropriate to assign to the neutralino local density p x the standard value 
for the total dark matter density pi = 0.3 GeV cm -3 , unless one specifically verifies 
that the neutralino relic abundance Cl x h 2 turns out to be at the level of an (fi/i 2 )min 
consistent with pi. This is why a correct evaluation of the event rates for % detection 
also requires a calculation of its relic abundance. 

Thus we evaluate Cl x h 2 and we determine p x by adopting a standard proce- 
dure [4]: when Q, x h 2 > (fi/2. 2 ) m in 5 we put p x = pi; when Q, x h 2 turns out to be less 
than (fi/i 2 )min5 we take 

fl x h 2 

Here (fi/i 2 )min is set equal to 0.03. 

For the neutralino relic abundance Q, x h 2 we employ the results of Ref.5. 
In Fig. 2 we report Cl x h 2 as a function of m x for three representative neutralino 
compositions: i) a gaugino-dominated composition (P = 0.9), ii) a composition of 
maximal gaugino-higgsino mixing (P = 0.5), hi) a higgsino-dominated composition 
(P = 0.1). As expected, out of the three compositions displayed in Fig. 2 the 
gaugino-dominated state provides the largest values of Cl x h 2 . In order to have more 
substantial values of Cl x h 2 , one has to consider more pure gaugino compositions 
(P > 0.99). This is explicitly shown in Fig. 3, where together with the reference value 

P = 0.5 also displayed are the rather extreme cases: P = 0.01 (very pure higgsino 
composition) and P = 0.99 (very pure gaugino composition). The very pronounced 
dips in the plots of Fig.s 2-3 are due to the s-poles in the X ~ X annihilation cross 
section (exchange of the Z and of the Higgs neutral scalars). 



(3) 



In Fig. 4 we display the values of Q, x h 2 versus m x in the form of a scatter 
plot obtained by varying M2 and fj, over a grid of constant spacing in the log-log 
plane of Fig.l. 

By comparing Fig. 4 with Fig. 3 one sees that, at fixed ra x , the minimum 
of Cl x h 2 is provided by the %-configuration of maximal mixing, as is expected; 
furthermore, one sees that the large spread in values for Cl x h 2 , displayed in Fig. 4, 
is due to configurations of extremely pure composition. 

3 Looking down, watching at muons that come up 



Let us turn now to the indirect search for neutralino dark matter which can be 
performed by means of neutrino telescopes [8] . Neutralinos, if present in our galactic 
halo as dark matter components, would be slowed down by elastic scattering off the 
nuclei of the celestial bodies (Sun and Earth) and then gravitationally trapped 
inside them. Due to the process of neutralino capture these macroscopic bodies 
could accumulate neutralinos which would subsequently annihilate in pairs. An 
important outcome of this x~X annihilation would be a steady flux of neutrinos 
from these celestial bodies. 

The differential neutrino flux at a distance d from the annihilation region is 
given by 

dN v _ T A ^ ( F) dN fv 
dE v ~ And 2 jr* xf dE v 

where Y a is the annihilation rate and F denotes the x~X annihilation final states 
which are: 1) fermion-antifermion pairs , 2) pairs of neutral and charged Higgs 

bosons, 3) one gauge boson-one Higgs boson pairs, 4) pairs of gauge bosons; 

denotes the branching ratio into the fermion / (heavy quark or r lepton), in the 
channel F; dNf v /dE v denotes the differential distribution of the neutrinos generated 
by the semileptonic decays of the fermion /. The z^'s, crossing the Earth, would 
convert into muons and generate a signal of up-going muons inside a neutrino 
telescope. Calculations of this muon flux from the original neutrino flux may be 
performed using standard procedures [8]. 

Particular care has to be taken in the evaluation of the annihilation rate Ta- 
This quantity is given by [9] 

T A = |tanh 2 ( — ) (5) 
Z Ta 

where t is the age of the macroscopic body (t = 4.5 Gyr for Sun and Earth), 
ta = (CCU) -1 / 2 , C is the capture rate of neutralinos in the macroscopic body and 
C a is the annihilation rate per effective volume of the body. The capture rate C is 
provided by the formula [10] 



(4) 



C= Px^_^L {MBf . ){v 2 te) . Xii (6) 



where p x and v x are the neutralino local density and mean velocity, <r e i,i is the cross 
section of the neutralino elastic scattering off the nucleus i of mass mi (for some 
properties of the elastic %-nucleus cross section see next Sect. 4 and Ref. [11]), 

Msfi is the total mass of the element i in the body of mass -Mb, { v 1 S c)i 1S ^ e 
square escape velocity averaged over the distribution of the element z, X{ is a factor 
which takes account of kinematical properties occurring in the neutralino-nucleus 
interactions. 



C A is given by [9] 



C A = <aV> ( m * ) 3 / 2 (7) 
Vo 40GeV ; K) 



where a is the neutralino-neutralino annihilation cross section and v is the relative 
velocity. V is defined as V = {3m 2 Pl T/(2p x 10 GeV)) 3 / 2 where T and p are the 
central temperature and the central density of the celestial body. For the Earth 
(T = 6000 K, p = 13 g • cm" 3 ) V = 2.3 x 10 25 cm 3 , for the Sun (T = 1.4 x 10 7 K, 
p = 150 g • cm" 3 ) Vo = 6.6 x 10 28 cm 3 . 

For the computation of the capture rate (and then also of ta) one has to 
use a specific value for the neutralino density p x . Here, for any point of the model 
parameter space, we evaluate p x as explained in Sect. 2. 

In Fig. 5 we give the results of our calculations for tanh 2 (^/r^) and for Ta, 
which include the rescaling for p x . This figure refers to the Earth for the case of the 
usual three representative %-compositions. For simplicity, in this figure, as well as 
in the following ones, only the results concerning positive values of the parameter 
p are shown; similar results hold for negative values of p. One clearly sees that 
equilibrium is not reached for m x > raw, because of the substantial suppression 

introduced in Ya by the factor tanh 2 (t / ta) ■ In Fig. 6 we show tank 2 (t / ta) and for 
Ta for the Sun; here the equilibrium between capture and annihilation is reached 
over the whole m x range. 

Now we report some of our results about the flux of the up-going muons 
in the case of x~X annihilation in the Earth. In Fig.s 7-8 we show the fluxes of 
the up-going muons as functions of m x for a number of values of the neutralino 
composition P, for x~X annihilation in the Earth. The threshold for the muon 
energy is = 2 GeV. We recall here that the present experimental upper bound 

for signals coming from the Earth is 4.0 • lO-^cm^s- 1 (90 % C.L.) [12]. By 
comparing this upper limit with our fluxes we see that the regions explored by 
Kamiokande (at our representative point: tan (3 = 8,mft, = 50 GeV) concern the 
mass range 50 GeV < m x < 65 GeV. These regions are illustrated in Fig. 9 in a 
p-M.2 plot. In this figure we also display the regions which could be explored by a 
neutrino telescope with an improvement factor of 10 (and of 100) in sensitivity. 



The location and the shape of the most easily explorable regions in the 



M.2 — fJ> plane depend on the Earth chemical composition and on the neutralino 
composition in terms of the gaugino, higgsino components. In fact the capture rate 
of neutralinos is more effective when neutralino mass matches the mass of some 
of the main components of the Earth and when the neutralino is a large gaugino- 
higgsino mixture. Because of these two properties the signal is maximal along iso- 
mass lines in the range 50-65 GeV, with elongations along iso-composition lines of 
sizeable mixing. 

In Fig. 10 we report the fluxes for up-going muons as functions of m x due 
to x~X annihilation in the Sun. As before the threshold for the muon energy is 
= 2 GeV. The evaluated fluxes are below the present experimental upper limit 
of Kamiokande: 6.6 • lO^cm"^" 1 (90 % C.L.) [12]. 

From these results it can be concluded that neutrino telescopes with an area 
above 10 s m 2 are very powerful tools for investigating neutralino dark matter in 
large regions of the parameter space. It also emerges from the previous results that 
the signals from the Earth and from the Sun somewhat complement each other to 
allow an exploration about DM neutralino over a wide range of m x . Taking into ac- 
count the appropriate on-source duty factor for the Sun, it turns out that the signal 
from the Sun starts overcoming the one from the Earth at about m x ~ 200 GeV. 
This occurs since, even if the Sun is mainly composed of very light elements, its 
gravitational field is much more effective compared to the Earth's one in capturing 
neutralinos. 



4 Patiently waiting that relic neutralinos hit our detector 



Another way to search for dark matter neutralinos is the direct detection 
which relies on the measurement of the recoil energy of nuclei of a detector, due 
to elastic scattering of %'s. The relevant quantities to calculate are the differential 
rate (in the nuclear recoil energy E r ): 

B- = N T ^r m " dvf{v)v d -^{v,E r ) (8) 
dE r m x Jv min (E r ) dE r 

and the integrated rate -Ri n t> which is the integral of Eq. (8) from the threshold 
energy which is a characteristic feature of the detector, up to a maximal energy 
£™ ax . In Eq.(8) Nt denotes the number of target nuclei, da e \/dE r is the differential 
elastic cross section and f(v) is the distribution of % velocities in the Galaxy. It 
is important to note again that the local density p x is evaluated here according to 
the procedure discussed in Sect. 2. In general, the %-nucleus cross section has two 
contributions: a coherent contribution, depending on A 2 (A is the mass number 
of the nucleus) which is due to Higgs and q exchange diagrams; a spin-dependent 
contribution, arising from Z and q exchange, proportional to A 2 J(J + 1). 

By way of example, let us remind the expression of the coherent cross section 
due to the Higgs-exchange [14]: 
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Figure 10. Fluxes of the up-going muons as functions of m x for x~X annihilation in the 
Sun, for the three representative neutralino compositions P = 0.1 (dotted line), 0.5 (short-dashed 
line), 0.9 (long-dashed line). The threshold for the muon energy is = 2 GeV. 
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where an is a quantity depending on the neutralino-Higgs and the neutralino- 
quarks couplings. It is worth mentioning that an depends rather sensitively on the 
%-composition and on a number of parameters, such as tan (3 and the Higgs masses. 

Except for very special points in the parameter space, the coherent contribu- 
tion to elastic cross section strongly dominates over the spin-dependent one. For a 
detailed analysis on the calculation of the direct event rates see [15] and references 
quoted therein. For an experimental overview about dark matter detectors see Ref. 
[16]. 

Here, as an example, we simply report in Fig. 11 the event rates i?i nt 
for a Germanium detector as a function of m x for neutralino compositions P = 
0.1,0.5,0.9. Rates are calculated by integrating the differential rate of Eq.(8) over 
the electron-equivalent energy range (2-4) KeV. Fig. 12 shows i?i nt for more pure 
neutralino compositions. In Fig. 13 we show the regions of the M2 — ft param- 
eter space which can be explored with an improvement of one and two orders of 



magnitude in the sensitivity of the detector. 

As for the shape of these regions we refer to the comments presented above, 
in Sect. 3, in connection with Fig. 9. Again, the signals are higher along the iso- 
mass line with an m x matching the mass of any element composing the detector. 
Thus, using detectors of different compositions allows explorations of the M2 — ft 
plane over a wide range in m x . For instance, investigation of regions with small m x 
values (of order of 10 GeV) with very low threshold detectors [17] would be very 
interesting. In fact this m x range (which is excluded by accelerator data only under 
a number of assumptions) is out of reach for the indirect detection discussed in the 
previous Section. 

5 An epilogue: how to save stationery and be happy 



The procedure for getting reliable results about neutralino relic abundance 
and detection event rates requires rather elaborate calculations, demanding great 
care and accuracy. This is because in sizeable regions of the parameter space many 
different channels and final states are competing, and the coupling constants of the 
various processes are very sensitive to the values of the free parameters. Furthermore 
interplay between calculations of different quantities has to be taken into account, 
if one really cares about realistic evaluations; for instance, this is the case when a 
rescaled neutralino local density has to be used in the evaluation of the event rates. 

It is clear that under simplifying assumptions and in special circumstances, 
for example when clear dominance of a particular process accidentally occurs, order- 
of-magnitude calculations may give reasonably good estimates. When this is the 
case, one may go through the straightforward exercise of feeding a few numbers in 
simple formulae (as the one in Eq.(9)) cheaply available in the market, and work 
out a reasonable (order of magnitude) result. It is obvious (but nevertheless worth 
mentioning) that these back-of-t he-envelope estimates, very enjoyable for pedagog- 
ical purposes, may in no way be interpreted as substitute for realistic calculations. 
The inadequacy of back-of-t he-envelope evaluations in this field may be proved 
quite easily, even on the back of a stamp if saving stationery and not caring about 
accuracy in calculations is the main point of the game! 
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